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Abstract 

In the iterative algorithm recently proposed by Waxman for solving eigenvalue problems, we 
point out that the convergence rate may be improved. For many non-singular symmetric potentials 
which vanish asymptotically, a simple analytical relationship between the coupling constant of the 
potential and the ground state eigenvalue is obtained which can be used to make the algorithm 
more efficient. 
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Recently, WaxmanQ] has proposed a convergent iterative algorithm for obtaining solu- 
tions of the eigenvalue problem, which does not involve a matrix diagonalization. For opera- 
tors which possess a continuum as well as a set of bound states this is a most advantageous Q] • 
In the case of ground state, for example, the eigenenergy, e, is determined numerically as a 
function of the coupling constant of the potential, A, and inverted to yield the e correspond- 
ing to the required value of A. The convergence rate of the algorithm, therefore, depends 
on two factors: the number of iterations required to find an eigensolution for a particular 
choice of e; and the number of times this must be repeated in order to determine the value 
of e which corresponds to the desired value of A. In this regard we wish to point out that for 
certain potentials a simple relationhip between A and e exists which may used to improve 
the convergence rate. 

In the Waxman algorithm , eigenpairs are determined as functions of the strength of 
the potential in the following manner. For simplicity consider a one- dimensional eigenvalue 
equation jl| 

\-d 2 x - XV{x)]u{x) = -eu(x) (1) 
lim u(x) = (2) 

|x| — >oo 

where d x = j-; A > is the strength parameter of the attractive potential (AV(x) > 
and V(x)— > as \x\ — > oo) and the energy eigenvalue, — e (with e > 0), is negative and 
corresponds to a bound state. Using Green's method a solution to eqQ is given by 

/oo 
G e (x-x')V{x')u(x')dx' (3) 
-oo 

where the Green's function G e (x) satisfies 

[-dl + e]G e (x) = 5(x) (4) 

lim G e (x) = 0. (5) 

\x\ — >oo 

Normalizing u(x) at an arbitrary x re f 

u(x ref ) = 1 (6) 

allows A to be written as (see eq(J3|)) 

A ~ J G e (x ref - x')V(x')u(x')dx' ^ 



which can then be used to eliminate A from eq(|3*|) 

f°° GJx - x')V(x')u(x')dx' 
u( x ) = ±22 e 2 1 y 1 y } (a) 

KJ fG e (x ref -x')V(x>)u(x')dx r {} 

Using equations (|7J) and (JHJ), A can be determined as a function of e in the following manner. 
For a particular choice of e eqtjHJ) can be iterated 

n J" G e (x re f — x')V (x')u n (x')dx' 

until it converges and A can then be determined from eq((Z|). Repeating for different values 
of e yields a set of different values of the potential strength A. When enough points have 
been determined, a simple interpolation procedure can be used to determine the dependence 
of e on A. 

On the other, a simple relationship between A and e can be obtained for non-singular 
symmetric potentials which vanish asymptotically in the following manner. Note that such 
the eigensolutions of such potentials have good parity. Taking the limit of eq(jT| 
yields 

linx^ (-dlu(x)) + e 

A = < 10) 

since u(x re f) = 1. If the potential is deep enough the ground state eigenvalue, e, is not 
small. In this case it is reasonable to let 

u( x ) = f( x )e- vrex x>0 (11) 

where a is constant and f(x) satisfies 

u(x) 0. (12) 

x— >oo 

The only dependence of u(x) on e is in the exponential tail of the wave function and is 
neglected in f(x) which generally has a polynomial structure. For small values of e neglecting 
this dependence is in most cases not possible. Then 

Jim. (-d 2 x u(x)) -»■ fixrefy-^**™* -2 y /alf'{x ref )e- VEEx ™f + ae (13) 
-> aie - ^"' + avy/ee-^"* + a 3 e (14) 

where ]' = %.■ Combining eqs((TDJ) and (fTl|) yields 

A = aie-^™' + a 2 Vee~ V ~ eXref + a' 3 e (15) 
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FIG. 1: The coupling constant A as a function of the ground state energy e for the lowest-lying 
even parity eigensolution of the inverted Gaussian potential (solid curve). The coupling constant 
for the potential is 1. The dashed curve is a fit to the data given by A = .21972 + 1.6087e 

Now for x re f = and potentials which are symmetric about 0, /' (0) = (=>- ci2 = 0) 
for the even parity solutions. In this case one obtains a simple linear dependence of A on 
e involving two unknown coefficients. In the more general case or for a different choice of 
x re f, three values of A must be calculated numerically for three different choices of e to 
determine the three coefficients in eq (|T3j) and the dependence of A on e. Care must be taken 
that eigensolutions used to determine A are well converged. Numerical errors may lead an 
incorrect dependence and require additional numerical determinations of A. 

In order to demonstrate the linear dependence of A on e, we have performed a calculation 
of lowest-lying even parity eigensolution of the inverse Gaussian potential|2| 

V(x) = e 2 . 

Note that x re f = and A=l. For each choice of e, A has been determined from eq(J7J) after the 
eigensolution, u(x), has converged. In figure 1 the linear dependence is clearly demonstrated 



and confirms numerically the existence of a simple linear relationship between A and e as 
long as e is not small. 
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